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Abstract. We give a detailed account of the techniques to compute radiative corrections in baryon 
semileptonic decays developed over the years by Mexican collaborations. We explain how the 
method works by obtaining an expression for the Dalitz plot of semileptonic decays of polarized 
baryons including radiative corrections to order ff(aq/nM\), where q is the four-momentum 
transfer and M\ is the mass of the decaying baryon. From here we compute the totally integrated 
0^ . spin angular asymmetry coefficient of the emitted baryon and compare its value with other results. 
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O INTRODUCTION 

o : 

Right after Pauli formulated the neutrino postulate, Fermi introduced the field theoretical 
treatment of the process n — > p + e~ + v e in the early 1930's y]]. This is the first known 
example of a V-theory. In the ensuing years other forms of nuclear /3 decays were 
observed, which prompted Gamow and Teller to formulate an A-theory. The concept 
of a new class of interactions, the weak interactions, had just emerged. 

Both the V and A theories were fused into the V — A theory by Sudarshan and Mar- 
shak and also independently by Feynman and Gell-Mann in the late 1950's, motivated 
by the discovery of parity non-conservation in weak interactions by Lee and Yang (the- 
oretically) and Wu and Telegdi (experimentally). 

In this context the weak interactions are described by an effective Lagrangian 
■&eff(x) = —(Gf/V2)jI(x)J^(x) +h.c, where Gf is the Fermi coupling constant and 
the weak current Jx (x) has the V — A structure. Jx (jc) can be separated into weak 
leptonic j\ (x) and weak hadronic 7^ (x) currents, namely, Jx (x) = j\ (x) + J 1 ^ (x) , where 
the leptonic current can be written directly in terms of the lepton fields whereas the 
hadronic one can be decomposed into parts having definite flavor SU(3) transformation 
properties and can be written in terms of quark fields I2D. 

The Lagrangian Jzf e ff however, faces many problems and cannot be taken as a self- 
consistent quantum field theory of weak interactions. Among other aspects, i) it is not 
renormalizable; ii) at high energies it leads to a violation of unitarity, i.e., it brings in 
probability non-conservation; and Hi) it has no room for neutral currents. 

With the advent of gauge theories, the cornerstone of the theory of weak interactions 
became the SU(2)xU(l) Weinberg-Salam theory, which currently possesses a quite 
impressive experimental evidence [3]. Further work, both theoretical and experimental, 
has finally yielded to the standard model of elementary particles, which embodies our 



knowledge of the strong and electroweak interactions. 

The purpose of this paper is to briefly review the achievements of the past thirty years 
of theoretical activity in baryon semileptonic decays (BSD) from the Mexican perspec- 
tive. We give a detailed account of BSD, focusing on techniques to the calculation of 
radiative corrections to observables which have been the major contribution of local 
research groups. The paper opens with a historical account of the development of the 
theoretical approach and continues with an application to the evaluation of radiative cor- 
rections to the spin-asymmetry coefficient of the emitted baryon. Numerical results are 
discussed afterwards. The paper closes with a brief summary and conclusions. 

RADIATIVE CORRECTIONS IN BSD: THE EARLY YEARS 

In dealing with the radiative corrections (RC) to BSD it is very important to have an 
organized systematic approach to deal with the complications that accompany them. 
These RC depend on an ultraviolet cutoff, on strong interactions, and on details of 
the weak interactions, other than the effective V — A theory. In other words, they have 
a model-dependent part [0, They also depend on the charge assignments of the 
decaying and emitted baryons. Their final form depends on the observed kinematical and 
angular variables and on certain experimental conditions. Over the years, our approach 
to the calculation of these RC has been to advance results which can be established in as- 
much-as possible once and for all. This task is considerably biased by the experimental 
precision attained in given experiments and by the available phase space in each decay. 

A systematic study of the calculation of RC in BSD was initiated back in 1980 
^\ following an approach originally introduced by Sirlin [0] for the electron energy 
spectrum in neutron decay and later extended to the decay of polarized neutrons tAm- 
An expression for the electron energy spectrum including radiative corrections was 
obtained in Ref. [@], which was accurate enough to allow experimental analyses to be 
performed in a model-independent fashion, provided hard bremsstrahlung photons are 
experimentally discriminated. At that time, those results were not directly applicable 
to the experiments being performed, which had lower statistics and made no provision 
to discriminate against hard photons. Later, in the same spirit, the RC to the differential 
decay rate of polarized neutral and charged baryons were presented in Ref. [9]. From that 
result, obtaining the RC to the decay rates and angular coefficients was straightforward. 
This was the first attempt to obtain RC to integrated observables. 

A step further was taken in 1987, when in order to improve previous results on the 
RC to the electron energy spectrum of baryon /3 decay [9], a new theoretical expression 
was obtained, this time including all terms of order aq/izMi where q is the four- 
momentum transfer and M\ is the mass of the decaying baryon II 1011 . The result is 
suitable for model-independent analyses of very-high-statistics experiments, without 
any experimental constraint on the detection of hard inner-bremsstrahlung photons. In 
Ref. f| 1 0f| it was shown that the bremsstrahlung contribution to BSD can be computed in 
a model-independent way up to terms of first order in q, by u sing the Low theorem 

In the meantime, new high- statistics experiments in BSD |ll|] made Dalitz plot (DP) 
measurements feasible and the application of RC was necessary. Much of the work had 
already been advanced in early calculations for the energy spectrum of the charged lep- 



ton and for the decay distribution of polarized baryons [SfilsLQll], so the new task was 
to adapt the results for the DP. In this respect, an expression for the DP of semileptonic 
decays of charged and neutral baryons including RC to order a and neglecting terms 
of order aq was introduced in Ref. U3\- The virtual RC presented no new complica- 
tions. In contrast, the bremsstrahlung RC became rather involved. The approximation 
implemented in that work was that the real photons are not observed directly but indi- 
rectly when the energies of the final baryon and charged lepton are found not to satisfy 
the final three-body overall energy-momentum conservation. In consequence, a detailed 
kinematical analysis for determining the integration limits over the photon variables was 
mandatory. Besides, the infrared divergence in the bremsstrahlung had to be identified 
carefully along with the finite terms that accompany it. A proper choice of variables 
yielded the integrations feasible. No doubt this paper marked an important path toward 
new results. The immediate application was to the computation of RC to the DP to order 
aq/lzM\, both for charged [| 1 3M and neutral Ill4ll decaying baryons. 

In summary, from the early years of research we learned that following the analysis 
of Sirlin of the virtual RC in neutron beta decay JUl and armed with the theorem of Low 
for the bremsstrahlung RC ^ one can show that the model-dependent contributions 
to both corrections (introducing new form factors) contribute to orders (u/Tl)(q/M\) n 
with n = 2,3,..., while orders n = 0, 1 lead to model-independent final expressions, 
because their model dependence is absorbed into the already existing form factors. 
The RC to BSD obtained to these latter orders are then suitable for model-independent 
experimental analyses and are valid to an acceptable degree of precision in the near and 
intermediate future: they will be useful in BSD involving heavy quarks. 

RADIATIVE CORRECTIONS TO BSD: THE RECENT YEARS 

There are six different charge assignments in BSD, namely, A~ — ► B°(l~V[), A — > 
5+ (7~vLA+ -> B°(l+v l ), A -> B-(l + Vi), A++ -> 5 + (/ + V/), and A+ -> B ++ (rv f ). 
In Ref. [|15n it was shown that it is not necessary to calculate each case separately. The 
final results of the last four cases can be directly obtained from the final results of the 
first two cases. This saves a considerable amount of effort, since only the first two cases 
need be calculated in detail. Such detail requires the choice of specific kinematic situa- 
tions. It is the DP that is normally studied experimentally. However, energy-momentum 
conservation may allow to discriminate events where photons are emitted carrying away 
energy such that the BSD is placed outside the so called three-body region (TBR) of 
the DP of non-radiative BSD. The events with those photons belong to what we have 
referred to as the four-body region (FBR). In addition, when the initial baryon is polar- 
ized the angular correlations between that polarization and the emitted baryon and the 
emitted charged lepton involve different RC, and it is not possible to obtain the final 
results of one correlation from the final results of the other correlation. Back in the mid- 
dle 1990's the systematic study of RC to BSD was complete to both orders n = and 
n = 1 for unpolarized decaying baryons so it was time to tackle new problems by con- 
sidering the polarization of either the decaying or emitted baryons. After some work, 
the analysis was finished to order n = for polarized decaying baryons, covering the 
spin-final baryon momentum and spin-final charged lepton momentum angular corre- 



lations rTTeL [Tvh. Those results were indeed useful for obtaining theoretical expressions 
for the angular spin asymmetry coefficients of the emitted baryon and charged lepton, 
respectively. Within our approximations and after producing some numbers, our results 
agreed well with others already published |l8|1 . 

Nowadays our goal has been extended to cover both the spin-final baryon momentum 
and the spin-final charged lepton angular correlations to order n = 1 of both neutral 
and negatively charged decaying baryons, restricted to the TBR. The former problem 
has been already solved |[l9|l whereas the latter is in progress. A further analysis will 
take into account the FBR contribution to the RC |[20ll and also the polarization of the 
emitted baryon. This will be done in the near future. Let us mention that up to this order 
our results will be useful in high-precision experiments involving heavy quarks. 

In the next sections we will describe how to apply our methods to get the RC to the 
TBR of the spin-final baryon momentum angular correlation to order n = 1 . Basically 
we borrow some recent results of Ref. Ill9n and analyze the case of a negatively charged 
decaying baryon in order to illustrate the procedure. 



VIRTUAL RADIATIVE CORRECTIONS 



For definiteness, let us consider the BSD 



A -^B + l + Vi, 



(1) 



and let A denote a negatively charged baryon and B a neutral one, so that / represents 
a negatively charged lepton and V/ its accompanying antineutrino. The four momenta 
and masses of the particles involved in process CD will be denoted by p\ = (E\,p{), 
Pl = (£ , 2,P2) 5 I = (E,l), and p v = (Zi^py), and by Mi, M2, m, and m v , respectively 
Ill9n . Further notation and conventions can be found in those references. We organize 
our results to explicitly exhibit the angular correlation §1 ■ p2, where §1 is the spin of A. 
This choice of the kinematical variables yields the angular spin-asymmetry coefficient 
of the emitted baryon, denoted hereafter by 0Cg. 

The uncorrected transition amplitude Mq for process (OQ) is given by 



Gy 

M = -^[uB(P2)W fl (pi 1 p 2 )u A (pi)][ui(l)0^ l v v (p v )] 1 



(2) 



where ua, ub, u\, and v v are the Dirac spinors of the corresponding particles and 



W^puPl) = fM 2 )y»+f2{q 2 )°^ + h{q 2 )^ 



+ 



q v 



,2n 



Mi 



Mi 



7s- 



(3) 



Here = 7^(1 + 75), q = p\— Pi is the four- momentum transfer, and fi(q 2 ) and 
gi(q 2 ) are the conventional vector and axial-vector form factors, respectively, which are 
assumed to be real in this work. 



The observable effects of spin polarization are analyzed by the replacement 



ua{p\) -^^{s\)u A {p\) 



(4) 



where T.(s\) = (1 — >§/i)/2 is the spin projection operator, and the polarization four- 
vector s\ satisfies s\ • s\ = — 1 and si ■ p\ = 0. In the center-of-mass frame of A, s\ 
becomes the purely spatial unit vector §i which points along the spin direction. 

The virtual RC can be separated into a model-independent part M v which is finite 
and calculable and into a model-dependent one which contains the effects of the strong 
interactions and the intermediate vector boson IflGjLdll]. This model-dependence can 
be absorbed into Mo through the definition of effective form factors, hereafter referred 
to as f[ and g\. The decay amplitude with virtual radiative corrections My is given by 



My = M + M v , 



(5) 



where 



M, 



a 
2% 



(6) 



The model-independent functions 4> and 4> ; contain the terms to order ff(aq/lzM\) 



II 1311 and they reduce respectively to and ; of Ref. Ill2ll . in the limit of vanishing 
(Xq/llM\. The second term in Eq. © can also be found in this reference. 

At this point we can construct the DP with virtual RC by leaving E 2 and E as 
the relevant variables in the differential decay rate for process After making the 
replacement © in ©, squaring it, summing over final spin states, we have 



EIM, 



K\ 2 



M 



Mi 2 



(7) 



spins 



spins 



spins 



The spin-independent contribution to My in Eq. ©, denoted here by M' y , was ob- 
tained to order ^(a/n) in Ref. II 1211 and later improved to order @{aqj %M\) in 



Ref. [13], so we will borrow the latter result. We now focus here on the spin-dependent 

part My^ to this order of approximation along the lines of Ref. Il6ll . where only terms 
of order G{aj%) were retained. The differential decay rate can be written as 



_ dE 2 dEdQ. 2 d% 

dTy = — — 7= M 2 mm v 

{2ny 



2 V 

spins 



M'v\ 2 



M 



W.2 



spins 



dT'-dT { "\ (8) 



where dT' v corresponds to the differential decay rate with virtual RC of unpolarized 
baryons given by Eq. (9) of Ref. [ 1311 , except that we have chosen here, without loss of 
generality, a coordinate frame in the center-of-mass frame of A with the z axis along the 
three-momentum of the emitted baryon, whereas in Ref. [| 1 3M the z axis was chosen along 



(s) 

the three-momentum of the emitted charged lepton. Similarly dT v ' can be obtained by 
standard techniques. Thus the decay rate with virtual RC is 



dT v = dQ \a ; + - (B\ 4> + Bf{&) - §1 • p 2 



AZ + -(B' 2 <i> + B' 2 W) 
it 



(9) 



where dQ. is a phase space factor and A' Q , B\, and B" are given in Ref. lfl3ll whereas Aq 
can be found in Ref. lllql . They all depend on the kinematical variables and are quadratic 
functions of the form factors. Equation © is the differential decay rate with virtual RC 
to order &{aq/ltM{). It is model-independent and contains an infrared divergent term 
in 4> which at any rate will be canceled when the bremsstrahlung RC are added. 

When including terms of order 0[aqj %M\) the model dependence of the RC shows 
up. For the virtual part, one can handle this by introducing effective form factors, f[ and 
g' t , in the uncorrected amplitude Mo, in such a way that [Jg] 

oc oc 

f' l (q 2 lP+ -l)=f l (q 2 ) + -a l (p + -l) 1 g[(q 2 lP+ ■ I) = gl (q 2 ) + -bi(p+ ■ I), 

71 71 

fk{<l 2 ,P+- l )=fk(q 2 ) + -ak, g' k (q 2 ,p+-l)=g k (q 2 ) + -b k , (* = 2,3) 

i.e., f[ and g\ have a new dependence in the electron and emitted baryon energies other 
than the ones in the q 2 dependence of the original form factors; this can be seen through 
a\ and b\, which are functions of the product p + ■ I = (p\ + p 2 ) • I- For the remaining 
form factors, within our approximations, a k and b% (k = 2, 3) are constant. 



BREMSSTRAHLUNG RADIATIVE CORRECTIONS 

In this section we turn to the emission of real photons in the process 

A^B + i + vi + y, (10) 

where A, B, £, and V/ denote the same particles as in the virtual case and 7 represents a 
real photon with four-momentum k = (a), k) . 

The process (fTOl) itself is strictly speaking a four-body decay whose kinematically 
allowed region is the joined area A + B depicted in Fig. [Q The distinction between these 
two regions has important physical implications. Region A is delimited by 

Ef n <E 2 < Ef ax , m<E<E m (11) 

where E m = (M 2 — M| +tn 2 ) /2M\ whereas region B is delimited by 

M 2 <E 2 <Ef n , m<E<E b (12) 

with Ejj = [(Mi — M2) 2 + m 2 ]/2(Mi — Mi). Finding an event with energies E and E 2 
in region B demands the existence of a fourth particle which must carry away finite 
energy and momentum. In contrast, in region A this fourth particle may or may not do 
so. Consequently, region B is exclusively a four-body region whereas region A is both 
a three-body and a four-body region. We will refer to them as the four-body and three- 
body regions (FBR and TBR), respectively. Our analysis of the bremsstrahlung RC will 
consider process (flOl) restricted to the TBR. 

The starting point will be to obtain the amplitude of process (flOl) with the spin effects 
included, retaining terms of order G(aqj nM\) by following the theorem of Low 
The amplitude for process (flOl) . Mb, is given in Ref. ||13|l and will not be repeated here. 



E 2 




FIGURE 1. Kinematical region (A + B) as a function of E and E% of the four-body decay (fTOb - The 
region A corresponds to what is referred to as the three-body region in this work. 



The bremsstrahlung contribution to the DP is obtained from the differential decay rate 

\I 2 mm v d 3 p 2 d 3 l d 3 p v d 3 k 
{2nf E 2 ~E E v 2co 



M 2 mm v d 3 p 2 d 3 l d 3 p v d 3 k v 2 4 



spins, e 

where the sum extends over the spins of the final particles and the photon polarization. 
In analogy to the virtual RC, the substitution © in £ | Mg| 2 of Eq. (fT3l leads to 

spins, e spins, e spins, e 

and therefore the differential decay rate dTs is 

dT B = dT' B -dTf. (15) 

Except for minor changes, the quantity dT' B in Eq. (IT5T ) corresponds to the 
bremsstrahlung differential decay rate for unpolarized baryons given by Eq. (42) 

of Ref. [il3ll . As for the spin-dependent term, dT B \ one can find further details about it 
in Ref. [fl9ll . The result can be cast into the form 

dTf = -d£lh -p 2 [B' 2 I (E,E 2 ) +C A ], (16) 

where Iq{E,E 2 ) contains the infrared divergent term iflfjll which cancels the one of 
its virtual counterpart and Ca in infrared convergent. Ca is presented in two forms in 
Ref. H. The first one is given in terms of triple integrals over kinematical variables of 
the photon and the second one is fully analytic. 



The full bremsstrahlung differential decay rate dTg is now constructed by subtracting 
jrL^ from dT' R . This dTg is added to dTy to obtain 



dT = dCl 



a, 



Aq-AqSi -p 2 + -[0/- 0//Si -p 2 ; 



K 



(17) 



where the functions 6/ and 6jj can be found in Ref. 

El- 
Expression (fTTT ) is model-independent and strictly only well-defined in the TBR of 
the kinematically allowed region. Although this equation is composed of rather lengthy 
expressions, it has been organized in such a manner that is easy to use. 



SPIN-ASYMMETRY COEFFICIENT a B 



The DP (TT71 ) so organized allows the calculation of a B , which is defined as 



N+-N- 
'N++N- 



(18) 



Here N + [N~] denotes the number of emitted baryons with momenta in the forward 
[backward] hemisphere with respect to the polarization of the decaying baryon. Appro- 
priate integration of Eq. (flTT) leads to 



a B 



B 2 + (a/7t)A 2 
Bi + (a/n)A 1 



(19) 



where 



B 2 
Bi 



p „c-max 



Tp pmax 



n-mm 
C 2 



A'odE 2 dE, 
A f dE 2 dE, 



A 2 
Ai 



pinax 
c 2 



j? pmax 



C 2 



0ndE 2 dE, 
0/dE 2 dE, 



and Aq, Ou,A' , and 0/ are defined in Ref. [fl9|l . 

Equation (TT91 ) is a model-independent analytic expression for ag, including radiative 
corrections to order &(aq/ %M\ ) . The uncorrected asymmetry coefficient 0Cg is obtained 
by dropping the terms proportional to aj% from this equation. 



We can evaluate a s in order to make a comparison with previous works |i8|,LL6|]. We 
use definite values of the form factors in order to compare under the same quotations, 
but this does not mean that our calculation is compromised to any particular values of 
them. Therefore we use f\ = 1.27, g\ = 0.89, and f 2 = 1.20 for the decay A — > pev, 
fx = 1, g, = -0.34, and f 2 = -0.97 for the decay IT -> nev, and f x =Q, gl = 0.60, 
and f% = 1 .17 for the process E~ — > Aev. 

First, we can evaluate (Xb(E,E 2 ) [the same Eq. (fT9l) without performing the double 
integrals] in several points of the DR This evaluation is presented in Table Q] for the 
process E~ — > nev, where the first part corresponds to aB{E,E 2 ) to order 0(a/n) from 
Ref. OH, the second part is as{E,E 2 ) to order &{aqj %M\) from this work, and the last 
part was computed in Ref. 1180 and reproduced here for comparison. 



TABLE 1. Percentage 8a B (E,E 2 ) with RC over the TBR in IT -» nev decay 
(a) to order 0{a/n) of Ref. um\ (b) to order 0{aqjnM\) of this work; and (c) 
computed in Ref. IU8I1 . 

a (a) 



0.8067 
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0.0 


0.0 


0.0 


0.0 
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CA T 
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A A 
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(b) 
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(c) 
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0.15 
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Next, we integrate numerically Eq. (fT9l) to obtain ag. This evaluation is displayed in 
Table [2l where the entries of the second column correspond to dg, the following column 
is the percentage differences defined as das = OCb — and the last two columns are 



reserved to comparisons with Refs. I116l . ll8h . There is a very good overall agreement. 



TABLE 2. Values of as and comparison with other works. 



Decay 



ao 8 as (this work) 8 as Ref. [16] 8o,b Ref. [18] 



A — > pev 
£~ — > nev 
E~ -»■ Aev 



58.6 -0.09 -0.2 -0.1 

66.7 0.05 0.1 0.0 
7.2 0.08 



DISCUSSION 



We have presented a short review of the situation (past and present) of the achievements 
in the computation of RC in BSD developed by Mexican research groups. We have 
shown how to apply them to the particular case of the spin-asymmetry coefficient of the 
emitted baryon and compared with other approaches. 

We can claim that the advancement in this topic has been important over the past 
years so that our understanding on the subject is now clear. Our approach to compute 
RC can be used in model-independent analyses for charm-baryon semileptonic decays 
to a high degree of precision. Even for semileptonic decays of baryons containing two 
charm quarks, they provide a good first approximation. Of course the problem is still 
open. More work is needed but the approach can be applied straightforwardly. 
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1. For an interesting historical review see for example T.-R Cheng and L.-R Lee, Gauge theory of 
elementary particle physics (Oxford University Press, Oxford, 1984) and references therein. 

2. A. Garcia and R Kielanowski, The Beta Decay of Baryons, Lecture Notes in Physics Vol. 222 
(Springer- Verlag, Berlin, 1985). 

3. S. Eidelman et al. [Particle Data Group], Phys. Lett. B 592 1 (2004). 

4. A. Sirlin, Phys. Rev. 164, 1767 (1967). 

5. F. E. Low, Phys. Rev. 110, 974 (1958). 

6. A. Garcia and S. R. Juarez W., Phys. Rev. D 22, 1 132 (1980); 22, 2923(E) (1980). 

7. T. Shann, Nuovo Cimento 5A, 591 (1971). 

8. A. Garcia and M. Maya, Phys. Rev. D 17 1376 (1978). 

9. A. Garcia, Phys. Rev. D 25, 1348 (1982). 

10. S. R. Juarez W., A. Martinez and A. Garcia, Phys. Rev. D 35 232 (1987). 

11. M. Bourquin et al. Z. Phys. C 21, 1 (1983); S. Y. Hsueh et al. Phys. Rev. D 38 2056 (1988). 

12. D. M. Tun, S. R. Juarez W. and A. Garcia, Phys. Rev. D 40, 2967 (1989). 

13. D. M. Tun, S. R. Juarez W., and A. Garcia, Phys. Rev. D 44, 3589 (1991). 

14. A. Martinez, A. Garcia and D. M. Tun, Phys. Rev. D 47, 3984 (1993). 

15. A. Martinez, J. J. Torres, A. Garcia and R. Flores-Mendieta, Phys. Rev. D 66, 074014 (2002). 

16. R. Flores-Mendieta, A. Garcia, A. Martinez and J. J. Torres, Phys. Rev. D 55, 5702 (1997). 

17. A. Martinez, J. J. Torres, R. Flores-Mendieta and A. Garcia, Phys. Rev. D 63, 014025 (2001). 

18. F. Gluck and K. Toth, Phys. Rev. D 46, 2090(1992). 

19. M. Neri, A. Martinez, R. Flores-Mendieta, J. J. Torres and A. Garcia, Phys. Rev. D 72, 057503 
(2005); 72, 079901(E) (2005). 

20. S. R. Juarez W., Phys. Rev. D 53, 3746 (1996); 55 (1997) 2889. 

21. R. Flores-Mendieta, J. J. Torres, M. Neri, A. Martinez and A. Garcia, Phys. Rev. D 71, 034023 
(2005); 70, 093012 (2004). 



REFERENCES 



